RELATIVISTIC ANTIHYDROGEN PRODUCTION 
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Antihydrogen has recently been produced in collisions of antiprotons with ions. While passing 
through the Coulomb field of a nucleus an antiproton will create an electron-positron pair. In rare 
cases the positron is bound by the antiproton and an antihydrogen atom produced. We calculate 
the production of relativistic antihydrogen atoms by bound-free pair production. The cross section 
is calculated in the semiclassical approximation (SCA), or equivalently in the plane wave Born 
approximation (PWBA) using exact Dirac-Coulomb wave functions. We compare our calculations 
to the equivalent photon approximation (EPA). 



I. INTRODUCTION 

Antihydrogen, the simplest bound state of antimatter, has first been produced and detected in 1995 at CERN in 
the Low Energy Antiproton Ring (LEAR) Q . The synthesis of an antiproton and a positron has been done by passing 
relativistic antiprotons through a xenon target. Only a few electron-positron pairs are produced in these collisions. 
In rare cases the velocities of the outgoing positron and antiproton are sufficiently close together, so that the particles 
join. In lowest order the formalism for calculating the production of antihydrogen is the same as the one used to study 
bound-free pair production in ion-ion collisions. Bound-free pair production is an important process as it is one of 
the reactions limiting the luminosity of heavy ion beams in high energy colliders Q . The equivalent-photon approach 
(Weizsacker- Williams method) has been used for the calculations of antihydrogen in H and for electron capture by a 
heavy ions in It is of interest to perform exact calculations in the framework of SCA or, equivalently, PWBA 

theory and compare them with these calculations, especially at low energies of the colliding particle. The bound-free 
pair production by bremsstrahlung has been considered in B and was found to be negligible. A recent calculation 
of Baltz shows that contributions of higher order effects reduce the cross section given in lowest order perturbation 
theory by a small amount for 7 — ► 00 M. 

II. TOTAL CROSS SECTION 

The total SCA cross section for pair production with electron capture in a heavy ion collision is given by 
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Throughout the paper we will set h = c = 1. The charge numbers of the projectile and target are denoted by Zp 
and Zt', the velocity of the projectile in the target rest frame is given by (3. The momentum transfer from projectile 
to target is q whose absolute value is given by s = \q\. The component of q in the direction of the projectile is 
q z = jj — Ef ^ Ez ■ The total energy of the bound electron is Ef] the one of the positron in the continuum \Ej\. (Please 
note that Ei is negative.) The third component of the Dirac matrices is a^. In our calculation tpf (— ip™/) is the 
Dirac-Coulomb wave function of a if-shell electron (/t/ = — 1, Ef = m e \J 1 — C 2 with ( = aZr and magnetic quantum 
number mj). ipi (= '0™*) is the wave function of an electron with negative energy Ei in the continuum describing the 
positron g. 

Because of charge-conjugation invariance the same formalism is used to calculate the production of relativistic 
antihydrogen in the bound-free process: p+ Zp — ► H(ls) + Zp + e~ . 
Using current-conservation we can write dl]) as: 
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It should be mentioned that (|l|) and (g) are not exactly equal if the wave functions are not exact eigenfunctions of 
the Dirac Hamiltonian . The vector (3± is perpendicular to q and defined by 

&=*-(4)f. (3) 
The solutions of the Dirac equation in the Coulomb field 

V(r) = -^, ( = aZ T (4) 
r 



are given by §|L§ 
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The angular dependence is expressed by the spin-angular functions 
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The radial functions g K and f K for the bound state are 



K \ - I, 1=3 + ^sgn(K) . (7) 
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with the Bohr radius denoted by ao- For the continuum we use the radial functions 
\fE, K {r)J \E-mJ n i v ; 

- 4 (fer+« 1 i? 1 ( 7 ( /t ) + if?j 2 7 ( K ) + l, 2ifcr)l (9) 
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In (||) and (^) 7 (k), fc, ?7 and TV are given by 

l{K) = ^^e, k=y/E>-ml, n = —, N = r(27(K) + 1} 
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We rewrite the expression in parenthesis in the matrix element of (^|) in the spherical basis (eo,e±i) with e*o = f 
and use the expansion of the vector plane waves into electromagnetic multipoles. From the orthogonality of the 
spherical harmonics we get the incoherent sum over the multipoles. After some algebra one obtains a relatively simple 
expression for the cross section. 
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The radial integrals are given by 
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and are evaluated quickly by the method presented in pi] . 

III. RESULTS 

We focus our attention on the process p + Zp — ► H(ls) + Zp + e~. In Fig. [I] and Table | the total cross section 
for Zp = Zp = 1 is given as a function of the Lorentz 7-factor of the projectile in the target rest frame. The cross 
sections are calculated numerically with ( |l2| ) - ( |l6| ) . 

FIG. 1. The solid line shows the total cross section a^ ot f° r antihydrogen production with a positron bound in the if-shell 
(Zp — Zt = 1) as a function of the Lorentz 7-factor of the projectile in the target rest frame. The dotted line presents the 
part of (12) containing T±. The cross section scales with Zp. 

TABLE I. The total cross section for antihydrogen production (Zp = Zp = 1) is given for different 7's of the projectile in 
the target rest frame for the Fermilab experiment ]l^ |. The .ftT-shell capture cross section is denoted by H(ls) and the cross 
section for capture into all shells by H(all). The contribution of all higher shells is estimated to 20% of the K capture |q|l3|. 
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We compare our results with experimental data measured in the low energy region for La 57+ target jL4j. Our 
calculations for iGshell bound-free pair production of La 57+ with projectile energies 0.405 GeV/u (7 = 1.43), 0.956 
GeV/u (7 = 2.026), and 1.3 GeV/u (7 = 2.40) yields a K /Z 2 P = 9.93*10" 7 b, 9.12*10~ 6 b, and 1.78*10" 5 b, respectively. 
Adding contribution from higher shells, which are assumed to be about 20% |qjl3]|, we find good agreement with the 
experimental results. We also compare our results with those of Becker [^5). Our results agree well with the results 
given there for 1 GeV amu -1 (7 = 2.07) and 15 GeV amu -1 (7 = 16.1). Our results for Zp = Zt = 1 are 



3.92 * 10 pb and <r^ ot — 2.55pb, respectively. 



IV. CORRESPONDENCE TO EPA 

To compare (n2) with the photoproduction cross section we rewrite it as: 
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with u) = Ef — Ei and q\ — s 2 — q 2 ; q z — af ot as a function of cu is shown in Fig. ||[ 



11 + \ -ft 



2?± 



(17) 



FIG. 2. The cross section for antihydrogen production as a function of ui. The solid line shows the cross section for the 
Lorentz factor 7 = 3 of the projectile in the target rest frame. The dashed line shows the cross section for 7 = 6 and the dotted 
line for 7 = 10. 



The S-matrix element for photo-induced bound-free pair production (also known as photo-induced JT-shell capture) 
in Coulomb gauge is given by 
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Using the multipole expansion for e^e we find the expression for the total cross section: 
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where Tj_ is the same as in (|14|), with s = u>. Here we denote the four-momentum of the photon by q. If in ( |17| ) Ti is 
omitted, we can define the photo- induced cross section for 'virtual' photons (see (|] Eq. (1)) as 
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The expression (17) can now be written as <r£; t (<7 7 * ). At this stage we can introduce the equivalent photon approxi- 



mation (EPA) . In this approximation the q 2 -dependence of is neglected and (|20| ) is replaced by the cross section 
for real photons (|l9|). Now the integral over q± would diverge and we must introduce a suitable cutoff. Now we can 
write the cross section in the equivalent photon approximation as 
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Ti(cj,g 2 ) as a function of the Lorentz-invariant variable Q 2 = —q 2 and to is given in Fig. |[ The momentum q 2 ^ is 
related to q 2 by 
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FIG. 3. T±_ as a function of Q 2 = — q 2 for fixed tj-values. The solid line shows T± for cu — 2.1, the dashed line for cu — 3, the 
dotted line for u> = 5 and the dash-dotted line for u) — 10. 
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FIG. 4. Ti as a function of Q 2 for fixed a;-values. The solid line shows T; for lo — 2.1, the dashed line for lo — 3, the dotted 
line for ui = 5 and the dash-dotted line for lo = 10. 



From Fig. || one can extract the cutoff parameter for the EPA calculation. For the range of w-values contributing 
significantly to the total cross section (see Fig. ||) T± is almost constant up to Q 2 max ~ 4m 2 . Beyond this value it 
falls off rapidly. Therefore q\ max is given by q 2 lax = Q 2 nax ~ i^) 2 - I n Fig- |we also show Ti(to,q 2 ) as a function of 
Q 2 = —q 2 and u which has a similar behavior as Tj_. In Fig. |we compare EPA results with the SCA calculation 



FIG. 5. The solid line shows the full SCA calculation results. The dotted line shows EPA results with q 
For the photo-induced cross section in EPA calculation we used the expression (0). 
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In the high energy region EPA fits correctly the total cross section. For 7 >> 1, the part of the integral for the 
total cross section in ([jj]) containing T; tends to a constant (see Fig. ||). Therefore for high energies of the colliding 
particles a to t = &epa + & const- For high 7- values we are in agreement with the EPA results of Q^|> Due to the 
w-dependence of the cutoff q\ max — Am 2 — (^) 2 the EPA formula is valid down to 7 = 3 . But in the low energy 
region (7 < 3) EPA fails. Further for small 7's the contribution of TJ to the total cross section can no longer be 
neglected but has been ignored in the derivation of the EPA formula ( pTf ) . 

FIG. 6. The solid line presents the part of (12) containing Ti. 

An analytical expression for the photo cross section ( |l9| ) can be obtained from Sauter's formula for the photoelectric 
effect Jlq] by means of crossing symmetry: 
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with £j = and fcj = ^/e| — 1. It is an approximation for aZ^ << 1 and relativistic velocities of the unbound 
electron. We compare the exact cross section (|l9|) with the analytic result (^3|) . Very good agreement is expected for 
Zt = 1- This is shown in Fig. 0. With (|l9|) it is also possible to calculate er 7 for aZx ~ 1 where (|2^ ) overestimates 
the cross section (see jl 17|). 



FIG. 7. The solid line gives the results of the exact calculations ( |l9| ) and the dotted line the approximated analytical 
expression (G3) for Zp = Zt = 1. 



V. CONCLUSION 



We calculate bound-free pair production in the semiclassical approximation (SCA) using exact Dirac-Coulomb wave 
functions. We compare results for ion-ion collision with experimental and theoretical data of P,^| JT^JT5[ | and find good 
agreement with their results. 

Our calculations have the advantage remaining valid also when the condition Za << 1 is not fulfilled, since exact 
Dirac wave functions are used. We give expressions for the cross sections already integrated over the angular variables 
of the free electron (or positron). The radial form factor integrals are evaluated quickly by means of recurrence 
relations as given in Cm. Therefore we are able to sum over many partial waves and this allows the extension of our 
calculations up to high values of 7. 

We compare the exact results which those of the equivalent photon approximation (EPA). We show that good 

agreement can be obtained already starting at 7 > 3 if one uses an w-dependent cutoff of qj_ max = -^Qmax ~~ (^y) 2 

with Q 2 nax ~ Am 2 . This justifies the cutoff chosen in at high 7- values. At high energies the contribution of 
the longitudinal photon tends to a constant. Therefore the total cross section is found to be of the form a to t — 

CEP A + O const- 
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